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ABSTRACT 

Context. Fragmentation-branching ratios of electronically excited molecular species are of first importance for the modeling of gas 
phase interstellar chemistry. Despite experimental and theoretical efforts that have been done during the last two decades there is still 
a strong lack of detailed information on those quantities for many molecules such as C„, C„H or C3H2. 

Aims. Our aim is to provide astrochemical databases with more realistic branching ratios for C„ (n=2 to 10), C„H (n=2 to 4), and 
C3H2 molecules that are electronically excited either by dissociative recombination, photodissociation, or cosmic ray processes, when 
no detailed calculations or measurements exist in literature. 

Methods. High velocity collision in an inverse kinematics scheme was used to measure the complete fragmentation pattern of elec- 
tronically excited C„ (n=2 to 10), C„H (n=2 to 4), and C3FL molecules. Branching ratios of dissociation where deduced from those 
experiments. The full set of branching ratios was used as a new input in chemical models and branching ratio modification effects 
observed in astrochemical networks that describe the dense cold Taurus Molecular Cloud- 1 and the photon dominated Horse Head 
region. 

Results. The comparison between the branching ratios obtained in this work and other types of experiments showed a good agree- 
ment. It was interpreted as the signature of a statistical behavior of the fragmentation. The branching ratios we obtained lead to an 
increase of the C3 production together with a larger dispersion of the daughter fragments. The introduction of these new values in the 
photon dominated region model of the Horse Head nebula increases the abundance of C3 and C3H, but reduces the abundances of the 
larger C„ and hydrocarbons at a visual extinction A v smaller than 4. 

Conclusions. We recommend astrochemists to use these new branching ratios. The data published here have been added to the online 
database KIDA (Kinetic Database for Astrochemistry, http://kida.obs.u-bordeauxl.fr I. 

Key words. Astrochemistry - Molecular data - ISM: clouds - ISM: individual objects: Horse Head Nebula 
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j_J 1. Introduction discovery in th e early 70s, many stud i es hav e been devoted to 

d " those species (IVan Orden & SavkallvL 1 1 9981) . mostly on their 

Carbon clusters and highly unsaturated hydrocarbons are structural and spectroscopic properties. The spectroscopic stud- 

among the molecules most often observed in the Inter Stellar ies were pushed mainly by the dramatica n y increasing perfor- 

Medium (ISM) (http://astrochermsttyjietK They are presentm mances of the obserV ation capabilities, while for the structural 

diffuse clouds ([Hobbs & Campbell, 1982; MaieretaU I2Q0J studies it was also driven by the increasing use of carbon-based 

Galazutdinov et al., 2002; Roueff et al., 2002), in dense clouds materials 
dBell et all 119991; Fosse et all 120011; iDickens et all 120001) . in 



circumst ellar envelopes d Cernicharo et all 120021). protostar en- The origin of carbon clusters and highly unsaturated hydro- 
carbons and their abundances in the ISM objects is still a puz- 
zling question. It is thought that they may come from gas-phase 



velopes {Sakai et al., 2008), planetary nebula dCernicharo et al 



2001) and in photon dominated regions (PDR) (Tev ssier et al 



2004 ). They are al so present around evolved (carbon) stars synthesis as well as be released into the gas-phase of carbona- 

(Hinkleetall |1988|), in the Titan ionosp here dLavvas et all ce0 us solid materials that are present in the ISM. Therefore the 

|2008|) or in comet commae (|Helbert et alj |2005|). Since their detailed investigation of their origin and abundances requires 

homogeneous and heterogeneous chemical reactions to be ex- 

Send offprint requests to: e-mail: chabot@ipno . in2p3 . fr plored. Whatever the carbon reservoir, those molecules are pro- 
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cessed in the gas phase by neutral-neutral and ion-neutral colli- 
sions as well as cosmic rays (CR), ultra violet (UV) photons and 
dissociative recombination (DR) from collisions with the ther- 
mal electrons. These last three mechanisms lead to highly elec- 
tronically excited species that may undergo fragmentation. The 
detailed dynamics and the fragmentation channels have thus to 
be investigated to enable their incorporation into various astro- 
chemical models. 

In gas phase chemistry, chemical networks for carbon 
clusters and highly unsaturated hydrocarbons synthesis have 
been proposed for a large variety of astr o physical conditions 
(iMillar et all 120001: lJura & KrotoL [l990t Black & DalgarnoL 
1977 t jHerbst & Leund. 1 19891: iTurner etafl, EOOOt iRollig et all 
2007; Gr eenberglll976l) . Nevertheless, despite the strong theo- 
retical and experimental efforts in structure and spectroscopy, 
many of the reaction rates still remain uncertain and may conse- 
quently limit the confidence in chemical models. 

It is convenient to write the reaction rate coefficient for a 
reaction A+B — > C+D as the product of a total reaction rate co- 
efficient k, ot and a branching ratio factor (BR) 



k[A + B C+D] = k tot [A + B (AB)*] xBR[(AB)* C + D]. 

(1) 

The reliability of the reaction rate may be improved when the 
total rate and BR are independently predicted. Most of the reac- 
tion rates are sensitive to the surrounding medium (temperature 
and photon spectrum for instance). In Eq. (1) this dependency is 
implicitly included in the total rate, whereas BRs are assumed to 
be constant. This assumption of constant BRs may however be 
inappropriate for some reactions as for example, neutral-neutral 
or ion-neutral reactive collisions where energy barriers can be 
present. For electronic excitation mechanisms followed by dis- 
sociation that occur in the ISM, the intermediate complex AB is 
appropriately described by an ensemble of molecules prepared 
in various highly exited electronic states. It has been pointed out 
several times that this physical situatio n fits t he requirements of 
statistical theory concep t s very well dHerbstl[l978l:lLeger et all 
1989; iBettens & HerbsA 1 19951) . In the particular case of small 
species, non-statistical mechanisms need also to be considered 
because the density of states is generally not high enough, even 
if high-energy excitations are involved. Non-statistical behav- 
ior may also arise in rare photodissociation situations, where a 
dissociative st ate is resonant with a discrete line from a local 
photon source dvanDishoeckLll988l) . 

In the first part of this paper we will report on measured BRs 
after electronic excitations that take place in a high- velocity col- 
lision experiment. In the second part we will introduce through 
the comparison between the different processes that are responsi- 
ble for the electronic excitations in the ISM, the idea of statistical 
universal BR for cosmic -ray-induced dissociation, photodisso- 
ciation, and dissociative recombination processes. We will then 
propose to correct BRs of current online databases (OSU Q, 
UMIST Q) when they are not resulting from measurements or 
detailed calculations. In a last part we will observe the influence 
of the new branching ratios on two typical astrochemical models 
that are aimed to simulate the dense clouds and the Horse Head 
PDR. 



1 http://www.physics.ohio-state.edu/ eric/research.html 

2 http://www.udfa.net/ 



2. Fragmentation of C n =2-io> C n=2 -4H and C 3 H 2 
molecules that are electronically excited by 
high-velocity collision 

2.1. Experimental set-up 

The experimental set-u p has already been described in de- 
tail in previous papers (IChabot et all l2002t iTuna et all [2008) 
and only a brief description is given here. Tandem MP (15 
M Volts) accelerator at the Institut de Physique Nuclaire d'Orsay 
prod uced molecular ca tionic beams of several MeV of en- 
ergy dWaast et all [1996). After magnetic analysis and collima- 
tion, singly charged molecules collided at high velocity (few 
a.u.) with a single target atom. Half a meter downstream, a 
parallel-plates electrostatic deflector analyzed parents and frag- 
ments with respect to their q/m ratio. A set of silicon detectors 
intercepted all the trajectories in a dedicated chamber. It is be- 
cause of their high velocity that parent and fragments may be 
detected with silicon detectors. This type of detector enables 
the measurement of the kinetic energy of the particle, i.e., with 
high-energy molecular beams of constant velocity, the mass of 
the detected molecule. Moreover, silicon detectors are 100% effi- 
cient and, owing to the kinematics, small detector sizes can cover 
100% of the solid angle for fragments emission. All the detectors 
are operating in coincidence, event by event. To get branching ra- 
tios on neutral species, the grid technique dMorgan et all 11971b 
iLarssonl Q995) may be used when the number of channels is 
small. In this method, a grid of known transmission is placed in 
front of the detector. Branching ratios of dissociation are linked 
to the recorded mass spectra through a set of linear equations. In 
the present work, because high-energy molecular beams were 
used, the shape of the current signal from silicon detectors could 
be analyzed to reso lve a pile-up of several neutral fragments 
dChabot et all 120021) . However, this technique was insufficient 
to get the whole information for the hydrogenated species. The 
grid technique was then mixed with a signal shape analysis to fi- 
nally fully resolve the fr agmentation of ne utral and cationic C n , 
C n H, and C3H2 species dTuna et alll2008l) . 

2.2. High-velocity collision (HVC) processes 
2.2.1. The excitation processes 

During the fast (~ 10" 16 s) collision between a molecule X and 
an atom, charge transfer (HVC - CT) may occur: 

X + + He ^ X* + He + (2) 

Owing to the high initial velocity of the transferred electron, 
electro nically excited states as well as the ground state are pop- 
ulated (Chab ot et all l2006l) . Until now no calculations have 
been performed to predict the internal energy distribution asso- 
ciated with this process. Nevertheless, internal energy distribu- 
tions may be deduced for the present species from experimental 
multiplicity distribution, i.e. probabilities a ssocia ted with a given 
number of emitted fragments dTuna et all 12008). Some internal 
energy distributions are reported in Fig.Q] 

In addition excitation (HVC - E) may also occur during the 
collision. 

X + + He ^ X + * + He (3) 

Let us note in t his velocity regime t hat only electronic ex- 
citations take place ( Wo hrer et all 1 19981) . Calculations of inter- 
nal energy following HVC-E have be en performed through in- 
dependent atom and electron model dWohrer & Watsorl [1993) 
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10 20 30 

Internal energy (eV) 

Fig. 1. Internal energy distributions following HVC-charge 
transfer (see text). Thick solid line: C2H, dotted line: C3H, bro- 
ken line: C4H, thin solid line: C3H2. Peak energies and widths of 
distributions may be regarded with errors on the ord er of a few 
ev. De tails on the method of extraction are given in Tun a et all 
(200§. 
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Fig. 2. Internal energy distributions following HVC-excitation 
(see text). Solid line: C2H + , dotted line: C3H + , broken line: 
C4H + , dot-dashed line: C3H}. Peak energies and widths of dis- 
tributions may be regarded with errors on the ord er of a few 
ev. De tails on the method of extraction are given in iTuna et aTl 
d2008h . 
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Fig. 3. Theoretical fragmentation BRs as a function of internal 
excitation energy for a C7 carbon cluster. Ca lculations were per- 
forme d with the MMMC statistical theory (Dfa z-Tendero et aTL 
2006). 
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Fig. 4. C7 calculated BRs with various internal energy distri- 
bution. The BRs were obtained by convolution of theoretical 
curves of Fig. [3] with the gaussian energy distribution: p(E) = 
: exp(-4(^-^) 2 ). The BRs are normalized to the multiplic- 



ity probability. Solid lines: <x = 2 eV dotted lines: cr — 5 eV.a) BR 
for two fragments breakup: triangles up: C4/C3, squares: C5/C2, 
circle: C(,/C. b) BR for three fragments breakup: hexagons: 
C3/C3/C, circle: C3/C2/C2, triangle down: C4/C2/C and C5/C/C. 



using doubly differential probabilities with the energy and im- 
pact parameter calculat ed within the c lassical trajectory monte 
carlo (CTMC) theory dMavnardlll999l) . These calculated distri- 
butions agree very well with the energy distributions shown in 
Fig. [2 wich are deduced from experimental multiplicity distri- 
butions. Ionization is also ocuring in HVC. In the present exper- 
iments with cationic beams it leads to multiply c harged species 
which are generally also electronically excited (IChabot et all 
2010). We will not report on the fragmentation of these multiply 
charged species, because this is out of the scope of this paper. 



2.2.2. Fragmentation of electronically excited species by 
high-velocity collision 

Statistical hypothesis has been proposed and often applied to cal- 
culate the fragmentation of finite size systems. It stipulates that 
all accessible micro-states are equiprobably populated, and con- 
servations of energy and momentum define which these micro- 
states are. Their ensemble is forming a so-called phase space. 
Theoretical predictions have to rely on numeration of individual 
states or integration in the phase-space. Two physical situations 
are presumably close to the statistical behavior. In the first one, 
the system gets enough time to explore all the phase space what- 
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Fig. 5. Experimental branching ratios of neutral C3H2 produced 
in HVC-CT reaction at a fixed number of emitted fragments (Nf) 
from left to right (N f ) = 2, 3, 4. 



Table 1. Probability of dissociation of C3H2 produced in HVC- 
CT reaction as a function of number of the emitted fragments 
(Nf).Nf=l corresponds to non-dissociative HVC-CT. 



N f 


Proba (± err) 


1 


0.16 (0.05) 


2 


0.43 (0.07) 


3 


0.28 (0.08) 


4 


0.08 (0.03) 


5 


0.05 (0.01) 
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Fig. 6. Experimental branching ratios of C3H2 produced in 
HVC-E reaction at a fixed number of emitted fragments (Nf) 
from left to righ (N f ) = 2, 3, 4. 



ever the bottlenecks in this random exploration. In the second 
physical situation, the system has no time to explore the whole 
phase space, but is prepared in such a large set of initial states 
that sampling of the phase space is presumably achieved by the 
entrance channels. Fragmentation of molecules that are electron- 
ically excited by HVC processes are clearly related to the second 
scenario and therefore can generally be interpreted in the frame 
of a statistical theory. For example, a statistical theory of frag- 
mentation, the microcano nical metropolis monte carlo (MMMC) 
dGross & Hervieuxl 1 1995b . has been applied to the carbon clus- 
ters C n and found to agree very w ell with the experimental mea- 
surements dMartineteTallllOOS). 

For C n , C n H, and C3H2 the creation of each fragment costs 
roughly the same amount of energy (4-7 eV El)- 

Therefore the 

fragmentation BRs exhibit specific energy domains associated 
to a given multiplicity (see Fig. [3}. Normalizing the branching 
ratios to multiplicity probabilities somehow vanishes the role of 
the input energy on BRs. To illustrate this point, we present in 
Fig. [4] the calculated BRs for C7. They were obtained by con- 
volution of the theoretical curves of the Fig. [3] with Gaussian 
shapes centered at energies from to 15 eVThe variations are 
found to be small, i.e. within ± 10% . It shows that as long as 
phase space scrambling is achieved, normalized BRs are mostly 
insentive to the internal energy distribution, unlike the multiplic- 
ity distributions. 

2.3. Experimental results and analysis 

Results for C3H2 fragmentation (multiplicity distribution: 
Table[Q BR: Fig.© and C 3 H+ (multiplicity distribution: Tabled 
BR: Fig. |6| have been recently obtained. Detailed physical 



3 C g/H 2 formation energ y is only 3 eV, but there is a barrier close to 
6 eV fcxonori et all 12008 



Table 2. Probability of dissociation of (C3H2)* produced in 
HVC-E reaction as a function of number of the emitted frag- 
ments (Nf ). Non-dissociative excited species are not detected in 
the experiment. 



N f 


Proba (± err) 


2 


0.34 (0.02) 


3 


0.43 (0.03) 


4 


0.19 (0.08) 


5 


0.05 (0.01) 



discussion on this species will be done elsewhere. Note that 
the detachment of one or two hydrogens is the most favor- 
able channel. Fragme ntation results for C n and C n H species can 
be found elsewhere (|Martinet et all 120041 iDfaz-Tendero et ail 
2006; iTuna et all |2008|) .A subset of these results are shown in 
Figs.Qto [9] Briefly, the most favorable channels were found to 
be H production for C n H and C3 production for C„. These chan- 
nels are indeed the most exothermic one s. For C n , C3 is a magic 
numb er because of the shell closure dVan Orden & Savk ally, 
1998). 



3. Application of HVC fragmentation branching 
ratios in astrochemical reaction networks 

3.1. Statistical fragmentation relevance in the context of ISM 
chemistry 



Three physical processes are involved in the creation of 
transient molecular electronically excited species in ISM 

They are HVC with c osmic rays 



dDulev & WilliamsL Il984h 

and secondary electrons dSpitzer & Tomaskol 1968b. photoab- 
sorpti on in the various local radiation fields ( Dulev & Will iams. 
1984), and recombination be tween molecu lar ions and thermal 
electrons in the local plasma dHerbstl 12003). The two latter pro- 
cesses are generally dominating, except perhaps in dark clouds 
where CR-i nduced secondary electrons may have a signi ficant 
contribution dPrasad & Tarafdarlll983l:lGredel et allll987l) . 

Because of the energy and momentum conservation laws the 
internal energy, in the electronic recombination process, should 
be close to the neutral species ionization potential (IP). In hydro- 
carbon molecules, the IP is always higher than the dissociation 
energy (DE). Then, fragmentation will occur rapidly for small 
size molecules, as compared t o radiative deexcitation . Acco rding 
to the IP value s (8 to 12 eVlVan Orden & SavkailvMl998h and 
DE (4 to 7 eV lDfaz-Tendero et all 120061: ITuna et1illl2008l) . the 
channels that lead to two fragments are expected to be the most 
populated. 

Statistical fragmentation may be invoked in the dissociative 
recombination process because there are always many electron- 
ically excited states close to the IP in carbon and hydrocarbon 
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Fig. 7. Comparison between a two fragments breakup after 
dissociative recombination (DR) (blue diamonds) and HVC- 
CT (red hexag ons - this work). DR-BR are ad apted from 
lEhlerding et ail d2004l) fo r C 7 H: lAngelova et all d2004bl) for 



CrH. lAngelova et all (2004a) for C 4 H and iHeber etafl 



(2006;) for C4. The parenthesis on Y-axis labels mean that hy- 
drogen may be localized on either of the two fragments in DR 
experiments. 



molecules ( Ivan Hemert & van DishoecH |2008). As a result a 
rather efficient sampling of the phase space should be performed 
in the entrance channels, and dissociative recombination BRs 
should be governed by a statistical fragmentation as in HVC. 
Figure 7 shows a comparison between HVC-CT and DR results 
for C4, C2H, C3H, C4H, and C3H2. Some of those results 
(C3H and C4) are not BR but a summation of BRs. Indeed, a 
resolution of C-H bond breaking is not always achieved in DR 
experiments nearby storage ring facilities due to the limitation 
of the grid method (see experimental section). The DR data 
agree with HVC-CT data within ± 10% on average. This small 
discrepancy arises because the internal energy distributions are 
quite different between HVC-CT (large distributions) and DR 
(narrow distribution around IP). An exception is seen in C2H, 
where the agreement is poor. It may come from a non-statistical 
behavior of the fragmentation for a triatomic molecule this 
small. 
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Fig. 8. Comparison between a two fragments -breakup after pho- 
todissociation (black triangles) and HVC-CT (red hexagons - 
this work) . The BRs for photodissociation are adapted from 
IChoi et all d2000h measurements. 
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Fig. 9. Comparison between a two fragments-breakup after pho- 
todissociation (triangles) and HVC-E (hexagons - this work). 
The B Rs for photodissociation are adapted from iGeusic et al.l 
(1986) 



In the ISM, molecules are photo-dissociated and photo- 
ionized because of the absorption of UV photons. These UV 
photons are produced by nearby stars. The interstellar standard 
radiation field (ISRF) e mbedding the gas in the Galaxy has been 
determined by Draine dDrainel [l978h . In star-forming regions, 
the UV flux is the sum of the ISRF plus the emission spectra 
of the nearby stars. UV photons are absorbed in the continuum 
by dust and via discrete lines of the most important molecules. 
Photo-absorption rates thus depend on the optical depth. In the 
interior of dense clouds, UV photons are produced by the excita- 



tion in Lyman and Werner electro nic states of molecular h ydro- 
gen by CR-secondary electrons dPrasad & Tarafdailll983l) . 

Taking into account the strong density of electronically ex- 
cited states for the c arbon and hydrocarbon molecules in the 6- 
1 1 eV energy range dvan Hemert & van DishoecH l2008h . many 
electronic transitions are expected to occur, resulting in a wide 
distribution of prepared excited states. Then fragmentation, as 
for HVC and DR, is expected to be governed by statistical be- 
havior. Figure 8 presents fragmentation BR for C4, C5, and Ce 
carbon clusters obtained by photodissociation and those from 
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Table 3. New HVC-SUBR for CR reactions. Total rates and BR 
of OSU-01-2007 ( http:www.eric/o.html) are also reported. We 
do not report on Cio because this reaction is not included in 
OSU-01-2007. 



Reactants 


Total rate 


Products 


HVC- 


OSU-01- 




(10 3 s- 1 ) 




SUBR 


2007-BR 




(OSU-01- 










2007) 








C 2 H + cr 


5.0 


C 2 + H 


0.81 


1.00 






C+CH 


0.19 


0.00 


C 3 H + cr 


5.0 


C 3 +H 


0.65 


1.00 






C 2 H + C 


0.33 


0.00 






C 2 + CH 


0.02 


0.00 


C 4 + cr 


1.0 


c 3 + c 


0.77 


1.00 






c 2 + c 2 


0.23 


0.00 


C 3 H 2 + cr 


5.0 


C 3 H + H 


0.46 


1.00 






C 3 +H 2 


0.24 


0.00 






C 2 H 2 + C 


0.23 


0.00 






C 2 H + CH 


0.04 


0.00 






C 2 + CH 2 


0.03 


0.00 


C 4 H + cr 


5.00 


C 4 +H 


0.58 


1.00 






C 3 H + C 


0.26 


0.00 






C 3 +CH 


0.00 


0.00 






C 2 H + C 2 


0.16 


0.00 


C 5 + cr 


1.00 


c 4 + c 


0.13 


1.00 






c 3 +c 2 


0.87 


0.00 


C 6 + cr 


1.00 


c 5 +c 


0.09 


1.00 






c 4 + c 2 


0.11 


0.00 






c 3 +c 3 


0.80 


0.00 


C 7 + cr 


1.00 


Cg+C 


0.01 


1.00 






c 5 + c 2 


0.19 


0.00 






c 4 + c 3 


0.80 


0.00 


C 8 + cr 


1.00 


c 7 + c 


0.03 


1.00 






c 6 + c 2 


0.01 


0.00 






c 5 + c 3 


0.90 


0.00 






C4 + C4 


0.06 


0.00 


C 9 + cr 


1.00 


c 8 + c 


0.00 


1.00 






c 7 + c 2 


0.06 


0.00 






c 6 + c 3 


0.66 


0.00 






c 5 +c 4 


0.28 


0.00 



Table 4. New HVC-SUBR for DR reactions. Total rates and BR 
of OSU-01-2007 ( http:www.eric/o.html) are also reported. We 
do not report about DR on C 2 H + and C+ because OSU-01-2007 
used experimental measurements. 



Reactants 


Total rate 


Products 


HVC- 


OSU-01- 




(10~ 7 cm 3 




SUBR 


2007-BR 




s- 1 ) (OSU- 










01-2007) 








C 3 H + + e~ 


3.0 


C 3 +H 


0.65 


0.50 






C 2 H+C 


0.33 


0.50 






C 2 + CH 


0.02 


0.00 


C 3 H+ + e- 


3.6 


C 3 H + H 


0.46 


0.42 






C 3 +H 2 


0.23 


0.42 






C 2 Ho + C 


0.24 


0.08 






C 2 H + CH 


0.04 


0.00 






C 2 + CH 2 


0.04 


0.08 


C 4 H + + e~ 


3.00 


C 4 +H 


0.58 


0.40 






C 3 H + C 


0.26 


0.15 






C 3 + CH 


0.00 


0.15 






C 2 H + C 2 


0.16 


0.30 


q + e" 


3.00 


C 4 + C 


0.13 


0.50 






c 3 + c 2 


0.87 


0.50 


q+e" 


20.00 


c 5 + c 


0.09 


0.50 






c 4 + c 2 


0.11 


0.50 






c 3 + c 3 


0.80 


0.00 


q + e" 


23.00 


c 6 + c 


0.01 


0.43 




c 5 + c 2 


0.19 


0.43 






c 4 + c 3 


0.80 


0.14 


q + e" 


20.00 


c 7 + c 


0.03 


0.50 






c 6 + c 2 


0.01 


0.50 






c 5 + c 3 


0.90 


0.00 






c 4 + c 4 


0.06 


0.00 


q + e" 


20.00 


c 8 + c 


0.00 


0.50 






c 7 +c 2 


0.06 


0.50 






c 6 + c 3 


0.66 


0.00 






c 4 + c 4 


0.28 


0.00 


CJo + e" 


20.00 


Cg+C 


0.01 


0.50 




c 8 + c 2 


0.01 


0.50 






c 7 +c 3 


0.70 


0.00 






c 6 + c 4 


0.03 


0.00 






C 5 +C 5 


0.25 


0.00 



HVC-CT In the re ported photodissociation experimental results 
dChoi et all l2000h . the incident photon energy has been varied 
up to 6 eV with some contribution of multi-photon absorption. 
The horizontal bars on these measurements correspond to a vari- 
ation of the BRs with the different photon energies. One observes 
agreement within + 10% between photodissociation and HVC- 
CT BR. Note that these photodissociation ex periments were fa - 
vorably compared to statistical calculations dChoi et all |2000). 
Figure 9 presents the same comparison for cationic C^ (n= 5 to 
10) species. Photodissociation was produced by s i ngle and run- 
away multiphoton absorption (ISowa et all 1 1 99 lb iGeusic et all 
1986). Therefore, those results have to be considered with cau- 
tion. Inside the above restrictions, the BRs from photodissocia- 
tion and HVC-E processes still agree also quite well. 

In view of the above comparisons, a universality of statistical 
BRs may be proposed. In order to take into account the particu- 
larity of the different electronic processes, an error bar (1 cr) of 
+ 10% seems reasonable to apply. 

3.2. New BRs for astrochemistry databases 

We propose here, to use the present complete set of BR ob- 
tained with HVC as statistical universal BRs (HVC-SUBR) for 



C„, C„H, and C3H2 when they are missing or guessed in ISM 
databases. Table [3] presents CR HVC-SUBR for C,„ C„H, and 
C3H2 together with OSU-01-2007 CR-BR. We do not report the 
Umist data base values because the BR are almost the same as 
those of OSU-01-2007. When reactions were not considered in 
OSU-01-2007, we do however used Umist06 rates. Note that 
the most recent version of the osu database (USU.2009) has the 
same BR as OSU-01-2007 for these species. All CR-BR from 
OSU- 01-2007 result from the pioneering estimates (Leun g~et all 
[1981 . In most cases, they used zero level statistical behavior to 
predict BR, i.e. all the dissociation goes only to the most exother- 
mic channel. It is consistent with the HVC-SUBR obtained for 
the C„H, but not for the C„ species. Indeed, they assumed C„_i/C 
and C„_i/C2 to be the most exothermic channels, whereas it is 
always the C„_3/C3 channel. In general, HVC-SUBR lead to a 
much wider dispersion of the fragmentation channels. It is re- 
markable that C3 production by CR is strongly enhanced. 

Table H presents new DR HVC-SUBR for C+, C„H + , and 
C3H+ together with OSU-01-2007 BR and reaction rates. We do 
not report on C? H + and Ct because OSU-01-2007 used e xper 



imental DR-BR dEhlerding etail l2004t iHeber et all 120061). For 
C3H + , in spite of existing partial measurements ( Angelova et all 



Chabot et al.: Statistical Universal Branching Ratios 



7 



Table 5. New HVC-SUBR for photo dissociation reactions. 
Total rates and BR of OSU-01-2007 (http:www.eric/o.html) are 
also reported. 



Reactants 


Total rate 


Products 


HVC- 


OSU-01- 




(10- 9 s" 1 ) 




SUBR 


2007-BR 




(OSU-01- 










2007) 








C 2 H + hv 


1.0 


C 2 + H 


0.81 


1.00 






C + CH 


0.19 


0.00 


C 3 H + hv 


1.0 


C 3 +H 


0.65 


1.00 






C 2 H + C 


0.33 


0.00 






C 2 + CH 


0.02 


0.00 


C 4 +hv 


0.4 


c 3 + c 


0.77 


0.50 






c 2 + c 2 


0.23 


0.50 


C 3 H 2 + hv 


2.9 


C 3 H + H 


0.46 


0.66 






C 3 + H, 


0.23 


0.34 






C 2 H 2 + C 


0.24 


0.00 






C 2 H + CH 


0.04 


0.00 






C 2 + CH 2 


0.04 


0.00 


C 4 H + hv 


2.0 


C 4 + H 


0.58 


0.50 






C 3 H + C 


0.26 


0.00 






C 3 + CH 


0.00 


0.00 






C 2 H + C 2 


0.16 


0.50 


C 5 +hv 


0.01 


c 4 + c 


0.13 


0.00 






c 3 + c 2 


0.87 


1.00 


C 6 +hv 


1.0 


c 5 + c 


0.09 


1.00 






c 4 + c. 


0.11 


0.00 






c 3 + c 3 


0.80 


0.00 


C 7 +hv 


1.0 


c 6 + c 


0.01 


1.00 






c s + c. 


0.19 


0.00 






c 4 + c 3 


0.80 


0.00 


C 8 +hv 


1.0 


c 7 + c 


0.03 


1.00 






c 6 + c 2 


0.01 


0.00 






c 5 + c 3 


0.90 


0.00 






C4 4" C4 


0.06 


0.00 


C 9 +hv 


1.0 


c 8 + c 


0.00 


1.00 






c 7 + c, 


0.06 


0.00 






c 6 + c 3 


0.66 


0.00 






c 5 + c 4 


0.28 


0.00 


Cio + hv 


1.14 


c 9 + c 


0.01 


0.82 






c 8 + c 2 


0.01 


0.17 






c 7 + c 3 


0.70 


0.00 






Cg + C 4 


0.03 


0.01 






C 5 + c 5 


0.25 


0.00 



l2004bl) . the OSU-01-2007 database uses the same probabil- 
ity for the emission of C or H. Although [Leung et alJ |l984) 
assu med 100% of H emission, equiprobability was proposed 
by iMillar et alJ (Q~988), based on more detailed calculations. 
The HVC-SUBR are indeed between these two extreme sit- 
uatio ns. For CtH + , OSU-0 1-2007 uses the experimental DR- 
BR (Ang elova et al ., 2004a) for the breaking of the C-C bonds 
and, without any available information on the ratio between 
emission of C or CH, OSU-01-2007 took it equal to 1. The 
HVC-SUBR results show that actually CH emission is very un- 
likely. For C^H^ OSU- 01-2007 uses the experimental DR-BR 
dAngelova et all l2004bl) for the C-C bonds. For the missing ex- 
perimental information on the proportion between H and H2 
emission when the carbon skeleton stays intact and the propor- 
tion between C and CH when it is broken, OSU-01-2007 still 
assumed equal probability. As for C3H + , HVC-SUBR show that 
CH emission is unlikely. For the proportion of evaporation of 
H and H2, equal probability is not too far from HVC-SUBR. 
It is remarkable that the HVC-SUBR give a ratio H 2 /H (35%) 



in qualitative agreement with the recent experiment of neutral- 
neutral reaction C+C 2 H 2 -> C 3 Hj dLeonori et all 120081) . For the 
C* species, as for CR, note that HVC-SUBR predict an enhanced 
production of C3 driven by DR together with a more extended 
dispersion in the carbon cluster-mass daughters compared to 
OSU-01-2007. 

Table |5] presents a new photon HVC-SUBR for C„, C„H, 
and C 3 H 2 together with OSU-01-2007 BR and reaction rates. 
Note that the OSU database has been developed mainly for cold 
dark clouds. Therefore, references about BR for photodissoci- 
ation processes are scarce. OSU-01-2007 Photodissociation BR 
for C 2 H and C3H used zero level statistical picture: only H emis- 
sion is allowed. For C4H, it assumed equal probability between 
emission of C 2 and H. For all C„H species HVC-SUBR predict 
H emission to be by far the dominant channel. For C3H 2 , OSU- 
01-2007 allows H or H 2 emission. Emission of C should also be 
an open channel. For the C„ series, as for the two previous pro- 
cesses, OSU-01-2007 assumed exclusive C emission. Exception 
arises fr om Cio where RRKM st atistical calculations were per- 
formed (iBettens & Herbsttll995l) . Again HVC-SUBR predicts a 
stronger enhancement of the C3 cluster production compared to 
OSU-01-2007. 

4. Effects of the new branching ratios on chemical 
model predictions for dense clouds 

4.1. Dense cloud model 

In order to test the effect of the new branching ratios on chem- 
ical model predictions for dense clouds, we u sed the Nahoon 
chem ical model developed by V. Wakelam ( Wakelam et al., 
2004http://www.obs.u-bordeauxl .fr/amor/VWakelam). This 
model follows the time-dependent chemistry of gas-phase 
species at a fixed temperature and density. We used Nahoon 
for a single spatial point (0D). We considered "typical" dense 
cloud conditions: a temperature of 10 K, an H density of 
2 x 10 4 cirT 3 , a visual extinction (Av) of 10 and a cosmic -ray 
ionization rate of 1.3 x 10~ 17 s _1 . For the initial conditions, we 
started with all the elements but H in the a tomic form with th e 
low-metal elemental abundances from iGraedel et al.l ([1982). 
The standard network used for this analysis is OSU-01-2007 
(http://www.physics.ohio-state.edu/~eric/research.html), which 
contains 452 species and 4430 reactions. We updated the 
branching ratios listed in Tables [3] El and [5] for 1) cosmic -ray 
dissociations of C„ (n = 4 to 9), C„H (n = 2 to 4) and C 3 H 2 , 2) 
dissociative recombinations of C* (n = 5 to 10), C„H + (n = 3 
to 4) and C3H2, and 3) photo-dissociations of C„ (n = 4 to 10), 
C„H (n = 2 to 4) and C3H 2 . We then let the system evolve over 
10 8 yr, when it reaches steady-state. 

4.2. Results 

Figure [10] shows the ratio between carbon-bearing species 
abundances computed with the updated database and those 
computed with the standard network OSU-01-2007. The new 
branching ratios modify the species abundances at two different 
times. The first one is a very early stage, before 10 4 yr, which 
is irrelevant for dense cloud chemistry. The second is much 
later, between 10 6 and 10 7 yr, but then the effect of the new 
branching rati os is less important. A t the typical dense cloud 
age of 10 5 yr (Wakelam et all 120061) . the new branching ratios 
are unimportant. Higher ag es have however been f ound for 
TMC-1 with other models dTerzieva & Herbstl 119981) . so the 
period around 10 6 yr is not without interest. After 10 6 yr, the 



5. Effects of the new branching ratios on chemical 
model predictions for photon-dominated regions 



a) 




10 2 10 3 10 4 10 5 10 6 10 7 10 8 
Time (years) 



b) 




CgH 
C 3 H 2 



10 2 10 3 10 4 10 5 10 6 10 7 10 8 
Time (years) 

Fig. 10. Evolution with time of the ratios of abundances com- 
puted with the updated database and with the standard net- 
work (OSU-01-2007) for a dense clouds for a) C„ and b) C„H m 
species, (see text) 



maximum effect is obtained for the largest molecules. The C„H 
molecules are an exception because C7H and CsH are more 
influenced than C9H. All the C-bearing species abundances are 
decreased by the new BR by a maximum factor of two. The 
weak sensitivit y of TMC-1 chemistry to BR has a lready been 
pointed out dMillar et all 119881: iHerbst & Led 1 19971) 



5.1. PDR model 

We tested the influence of these new rates on a PDR model. For 
this, we u sed the Meudon PD R code (http//pdr.obspm.fr) de- 
scribed in iLe Petit et all d2006l) . The Meudon PDR code com- 
putes the structure of a ID plan-parallel and stationary slab of 
dust and gas. It consistently solves the radiative transfer from far 
UV to sub-millimeter, chemistry, and thermal balance. To test 
the influence of the new rates, we reproduce d the mode l of the 
Horse Head by this PDR code presented in iPetv e^ l. (2005) 
first with an old chemical network based on the rates provided 
by OSU-01-2007, secondly with the new branching ratios pre- 
sented in this paper. The Horse Head is good candidate because 
of the large number of observed hydrocarbons. Pety et al. (2005) 
suggest that fragmentation of PAHs can contribute to the synthe- 
sis of small hydrocarbons and conclude by mentioning the need 
for precise chemical rates to perform chemical models in these 
regions. The proton density in the Horse Head is estimated to 
njj = 10 s cm -3 and the intensity o f the incident UV flux to 100 
times the ISRF (in Draine's units, Draine, 1978). We adopted a 
flux of comic rays of 5 x 10~ 17 per H and per second. The model 
assumes a semi-infinite cloud. Our objective is to compare the 
profile of abundances of some relevant species computed with 
our new branching ratios and with older ones. To refine models 
of the Horse Head is beyond the scope of this paper. 

5.2. Results 

Figure [TT] presents the ratios of abundances provided by the 
two chemistries of C„ and C„H with n from 2 to 9 as a function 
of the position in the cloud expressed in visual extinction. 
This ratio is defined as the abundance provided with the new 
chemistry divided by the abundance obtained with a older 
rates. First we note that the effect of the new rates are only 
visible for Ay lower than 4. Indeed, it is in this region that the 
dissociative recombinations and photo-processes dominate the 
other chemical processes. 

It is often difficult to reproduce the abundance of C3 in 
PDR models (for example in the diffuse i nterste llar gas towards 
Zeta Persus, the model by ILe Petit et al.l d2004l) requires a high 
density component to reproduce this molecule). Figure [TT] 
shows that the new branching ratios enhance the abundance of 
this molecule. This is explained by two factors. First, the new 
branching ratios of the CsH + recombination reaction enhance 
the route leading to C3 by 30%. Secondly, we showed that 
the dissociative recombination of Cg can efficiently produce 
two C3 molecules. This route was not considered in previous 
chemistries. 

For C„ molecules with n >4, the new branching ratios systemat- 
ically produce a decrease in the abundances. The model shows 
that the abundance of C3H is significantly enhanced in the 
photodissociation area (Fig.fTTb). This is explained by the new 
photodissociation route of C4H. This route was not considered 
in the OSU database. 

Finally, the abundances of the hydrocarbons, C„H molecules 
(n>3), are reduced compared to the old chemistry. This is di- 
rectly linked to the decrease in the abundances the C„ molecules. 
Indeed, the chain of reactions leading to the hydrocarbons with 
n>3 is 
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Cn ^ C+ +1 * C n+I H + S C n H. (4) 

Because the abundance of C„ molecules is reduced by 
the new branching ratios, the abundances of the C„H are also 
reduced. 



6. Conclusions 

High velocity collision in inverse kinematics scheme was used 
to measure the complete fragmentation pattern of electronically 
excited C„ (n=2 to 10), C„H {n-2 to 4) and C3H2 molecules. 
Branching ratios of dissociation were deduced from those ex- 
periments. The comparison between the branching ratios ob- 
tained in this work and other types of experiments showed a 
good agreement. It was interpreted as the signature of a statisti- 
cal behavior of the fragmentation. We thus propose new branch- 
ing ratios for: 1) the dissociation of C„ (n=2-10), C„H (n-2-4) 
and C3H2 molecules by interstellar UV photons or cosmic-ray 
processes and 2) the dissociative recombination of (n=5-10), 
C„H + (n=3-4) and C3H2. The new values have been tested in 
dense cloud and PDR models. We showed that only chemistry 
occurring at Ay smaller than 4 is really affected. We however 
recommend astrochemists to use these branching ratios, even for 
dense chemistry, because it is well known that the importance of 
a reaction depends on the network we use. The data published in 
this paper have been added to the online database KIDA (Kinetic 
Database for Astrochemistry, http://kida.obs.u-bordeauxl.fr). 
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Fig. 11. Evolution, with the visual extinction Av, of the ratios of 
aboundances computed with the updated database and with the 
standard network (OSU-0 1-2007) for a PDR a) C„ and b) C„H m 
species (see text). 



